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Compelling evidence for a new form of matter has been claimed to be formed in Pb+Pb 
collisions at SPS. We critically review two suggested signatures for this new state of 
matter: First the suppression of the J/\l/, which should be strongly suppressed in the 
QGP by two different mechanisms, the color-screening []I| and the QCD-photoeffect 0. 
Secondly the measured particle, in particular strange hadronic, ratios might signal the 
^ . freeze-out from a quark- gluon phase. 

cn 

^ ■ 1. J/\l/ suppression 

o ■ 

O ■ The QCD factorization theorem is used to evaluate the PQCD cross sections of heavy 
quarkonium interactions with ordinary hadrons. However, the charmonium states (here 
denoted X) are not sufficiently small to ignore nonperturbative QCD physics. Thus, 
we evaluate the nonperturbative QCD contribution to the cross sections of charmonium- 
nucleon interaction by using an interpolation between known cross sections P]. The 



O 



C ■ J/'^-N cross section evaluated in this paper is in reasonable agreement with SLAG data 

i- 

I Indeed, the A-dependence of the J/\l/ production studied at SLAG at ~ 20 GeV 
I exhibits a significant absorption effect |Q leading to aabs{J 1^-^) = 3.5 ± 0.8 mb. It was 
demonstrated that, in the kinematic region at SLAG, the color coherence effects are 
still small on the internucleon scale for the formation of J/^&'s. So, in contrast to the 
findings at higher energies, at intermediate energies this process measures the genuine 
J/'^-N interaction cross section at energies of ~ 15-20 GeV 

To evaluate the nonperturbative QCD contribution we use an interpolation formula 
for the dependence of the cross section on the transverse size 6 of a quark-gluon config- 
uration Three reference points are used to fix our parametrization of the cross sections 
(cf. Tab. |l|). The X-N cross sections is calculated via: cr = /cr(6) ■ \^ {x,y , z)\'^dx dy dz, 
where '^{x,y,z) is the charmonium wave function. In our calculations we use the wave 
functions from a non-relativistic charmonium model (see [||). 

We follow the analysis of [0 to evaluate the fraction of J/\l/'s (in pp collisions) that 
come from the decays of the x and So, the suppression factor S of J/\I''s produced in 
the nuclear medium is calculated as: 

^ = 0.6 ■ (0.92 ■ 5-^/* + 0.08 ■ 5*') + 0.4 ■ S^. Here are the respective suppression factors 
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20.0 
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Table 1 

The total quarkonium-nucleon cross sections o. For the x two values arise, due to the spin 
dependent wave functions (Zm = 10, 11). 



of the different pure charmonium states X in nuclear matter. The S are for minimum 
bias J) A collisions within the semiclassical approximation (cf. [§). 

The charmonium states are produced as small configurations, then they evolve to their 
full size. Therefore, if the formation length of the charmonium states, //, becomes larger 
than the average internucleon distance, one has to take into account the evolution of the 
cross sections with the distance from the production point 0. 

The formation length of the J/\& is given by //• ~ ^ — , where p is the momentum 

of the J/\E' in the rest frame of the target. For a J/\E' produced at midrapidity at SPS 
energies, this yields If ^ 3 fm. Due to the lack of better knowledge, we use the same 

^ 3 fm for the x- For the \E'' we use If ^ Q fm, because it is not a small object, but has 
the size of a normal hadron, i.e. the pion. For Eiab = 800 AGeV we get a factor of two 
for the formation lengths due to the larger Lorentz factor. 

However, this has a large impact on the \E'' to J/\E'-ratio depicted in Fig. |T], which shows 
the ratio 0.019 ■ Sq,'/Sj/q, calculated with (squares (200 GeV) and triangles (800 GeV)) 
and without (crosses) expansion. The factor 0.019 is the measured value in pp collisions, 
because the experiments do not measure the calculated value Sxi,'/Sj/^ but ■ -^m/^ 

are the branching ratios for J/^E', ^1^' — > /i/x. 
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Figure 1. 

The ratio 0.019 ■ Sqii/Sj/^, is 
shown in pA (crosses) in compar- 
ison to the data (circles). The 
squares and the triangles shows 
the ratio calculated with the ex- 
pansion of small wave packages 
(see text). 



The calculations which take into account the expansion of small wave packages show 
better agreement with the data (circles) (taken from [^]) than the calculation without 
expansion time, i.e. with immediate J/'^ formation. If = {]. We calculated this effect 
both at Eiab = 200 AGeV and 800 AGeV. The data have been measured at different 
energies {Eiab = 200, 300, 400, 450, 800 GeV and ^/s = 63 GeV). One can see that this 
ratio is nearly constant in the kinematical region of the data, but it decreases at smaller 
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momentum (e.g. Eiab = 200 AGeV and y < 0) due to the larger cross section of the 

However, the P-states yield two vastly different cross sections (see Tab. ^ for xio and 
Xii, respectively. This leads to a higher absorption rate of the xu as compared to the 
Xio- This new form of color filtering is predicted also for the corresponding states of other 
hadrons; e.g. for the bottomium states which are proposed as contrast signals to the 
J/^'s at RHIC and LHC! 

Furthermore it is important to also take into account comoving mesons. Therefore 
we use the UrQMD model |10,11|. Particles produced by string fragmentation are not 
allowed to interact with other hadrons - in particular with a charmonium state - within 
their formation time (on average, Tp ^ 1 fm/c). However, leading hadrons are allowed to 
interact with a reduced cross section even within their formation time . The reduction 
factor is 1/2 for mesons which contain a leading constituent quark from an incident nucleon 
and 2/3 for baryons which contain a leading diquark. 

Figure |^ shows the J/ ip to Drell-Yan ratio as a function of Erp for Pb+Pb interactions at 
160 GeV compared to the NA50 data The normalization of 5^/,(t( J/^)/a(DY) = 

46 in pp interactions at 200 GeV has been fit to S+U data within a geometrical model 
1^. The application of this value to our analysis is not arbitrary: the model of Ref. ^ 
renders the identical i^y-integrated J/ijj survival probability, S = 0.49, as the UrQMD 
calculation for this system. An additional factor of 1.25 [T^ has been applied to the 
Pb+Pb calculation in order to account for the lower energy, 160 GeV, since the J/ip and 
Drell-Yan cross sections have different energy and isospin dependencies. 

The gross features of the Et dependence of the J/ip to Drell-Yan ratio are reasonably 
well described by the model calculation. No discontinuities in the shape of the ratio as 
a function of Et are predicted by the simulation. The new high Et data []T3| decreases 
stronger than the calculation. This could be caused by underestimated fluctuations of the 
multiplicity of secondaries in the UrQMD model. This occurs, since high i^^-values are a 



trigger for very central events with a secondary multiplicity larger than in average |15 




Figure 2. 

The ratio of J/ip to Drell-Yan 
production as a function of Et for 
Pb+Pb at 160 GeV (see text). 
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2. Particle production 

Ideal gas model calculations have been used for a long time to calculate particle produc- 
tion in relativistic heavy ion collisions, e.g. |T^-^. Fitting the particle ratios as obtained 
from those ideal gas calculations to the experimental measured ratios at SIS, AGS and 
SPS for different energies and different colliding systems yields a curve of chemical freeze- 
out in the T — /i plane. Now the question arises, how much the deduced temperature and 
chemical potentials depend on the model employed. Especially the influence of changing 
hadron masses and effective potentials should be investigated, as has been done for exam- 



ple in [22-25|. This is of special importance for the quest of a signal of the formation of 



a deconfined phase, i.e. the quark-gluon plasma. As deduced from lattice data the 
critical temperature for the onset of a deconfined phase coincides with that of a chirally 
restored phase. Chiral effective models of QCD therefore can be utilized to give important 
insights on signals from a quark-gluon plasma formed in heavy-ion collisions. 

Therefore we compare experimental measurements for Pb-|-Pb collisions at SPS with 
the ideal gas calculations and results obtained from a chiral SU(3) model [^,^]. This 



effective hadronic model predicts a chiral phase transition at T ^ 150MeV. Furthermore 
the model predicts changing hadronic masses and effective chemical potentials, due to 
strong scalar and vector fields in hot and dense hadronic matter, which are constrained 
by chiral symmetry from the QCD Lagrangean. 

In Ijl^ the ideal gas model was fitted to particle ratios measured in Pb+Pb collisions at 



SPS. The lowest is obtained for T = 168MeV and fiq = 88.67MeV. Using these values 
as input for the chiral model leads to dramatic changes. As can be seen in Figure ^ 
There are two main reasons for the strong deviations. First, since the chosen temperature 



O ig-PBM (T=168 MeV, /j„=88.67 MeV, ;j,=17.56 MeV) 
taken from Braun-Munzinger et al, PLB 465(1999)15 

A chiral (T=168 MeV, /i,=88.67 MeV, (i,= 17.56 MeV) 
both inch feeding 
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Figure 3. 
Comparison of 
culation with 
168MeV, /iq 



a 



chiral cal- 
for T = 
88.67MeV. The 



huge differences result from the 
dropping baryonic masses, de- 
pending on the strangeness con- 
tent of the particles, and from the 
change of the effective potentials 
in the chirally restored phase. 
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c; 



lies above the chiral phase transition temperature of the model, the effective masses of 
the baryons are lowered dramatically (see Figure ^ . The second reason are the strongly 
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Figure 4. 

Baryon masses as function of 
temperature for vanishing chem- 
ical potential One sees that a 
phase-transition occurs at Tc ~ 
150MeV. There the masses drop 
significantly and then nearly sat- 
urate. The baryon mass above 
scales with the strangeness con- 
tent of the corresponding baryon. 
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changed effective potentials. While the chemical potential of the proton with the chosen 
parameters is fip = 266MeV in the ideal gas model, the effective chemical potenial of 
the proton in the chiral model is only fip ^ 8MeV, due to the reduction by the vector- 
field. Such a small effective chemical potential leads to nearly equal particle-antiparticle 
numbers and therefore to tremendeously enhanced particle-antiparticle ratios. 

Obviously, the freeze-out temperature and chemical potential have to be readjusted to 
account for the in-medium effects of the hadrons in the chiral model. 

We call the best fit the parameter set that gives a minimum in the value of x^? with 



X 



[exp _ model\ ' 



(1) 



Here r^^^ is the experimental ratio, r^"'^'^'- is the ratio calculated in the model and cxj 
represents the error in the experimental data points as quoted in We included the 
following ratios in the fit procedure: ^, ^, f^,^, f, |, Ig, ^±^, The 



resulting values of for different T — fi pairs are shown in figure ^. In all calculations 
lis was chosen such that the overall net strangeness fg is zero. The best values for the 
parameters are T = 144MeV and fig ~ 95MeV. While the value of the chemical po- 
tential does not change much compared to the ideal gas calculation, the value of the 
temperature is lowered by more than 20 MeV. Furthermore Figure |^ shows, that the 
dropping effective masses and the reduction of the effective chemical potential make the 
reproduction of experimentally measured particle ratios as seen at CERN's SPS within 
this model impossible for T > T^. Only for temperatures below we obtain reasonable 
fits. Using the best fit parameters we compare the particle ratios, as calculated in the 
chiral model, with the ideal gas calculation of |18[ and data as compiled in |]18|. This is 
shown in figure ^ One obtains a satisfactory agreement over several magnitudes. The 
values of the chiral model is xlhirai — 26.5. This is larger than the value in the ideal 
gas model of [|n 



(xL = 13). Note that in ||T8[ weak decays are accounted for. If we 



use our ideal gas calculation without feeding from weak decays we obtain only slightly 



changed best values (T = 168, /Xg = S2MeV) compared to and xlg-PFM = 21.6 



This shows that the chiral and ideal gas analysis using the same feeding procedures yield 
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Figure 5. 

data taken 
The best fit parame- 
ters are T = lAAMeV and /ig ~ 
95MeV. fig is chosen such that 
fs = 0- For T > Tc no agreement 
with data from CERN SPS lead 
on lead collisions can be obtained. 
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comparable agreement with data concerning the value of x^- Figure ^ shows that there 
is satisfactorily agreement between data and experiment. We want to emphasize, that in 



O SPSPb+Pb 

taken from Braun-Murizinger et al, PLB 465(1999) 15 
A chiral T=144 MeV, ;i,=95 MeV 
□ ideal gas calculation, PLB 465(1999) 15 
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Figure 6. 

Comparison of chiral model for 
T = 144 MeV compared to data 
and ideal gas results from 
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spite of the strong assumption of thermal and chemical equilibrium the obtained values 
for T and fi differ significantly depending on the underlying model, i.e. whether and how 
effective masses and effective chemical potentials are accounted for. In our model, we 
conclude that the observed particle ratios as measured at CERNs SPS do not signal the 
freeze-out from a chirally restored quark-gluon phase. Note that we assume implicitly, 
that the particle ratios are determined by the medium effects and freeze out during the 
late stage expansion - no flavor changing collisions occur anymore, but the hadrons can 
take the necessary energy to get onto their mass shall by drawing energy from the fields. 
Rescattering effects will alter our conclusion but are presumably small when the chemical 
potentials are frozen. 
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